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ABSTRACT 

It is widely believed that T Tauri winds are driven magnetocentrifugally from 
accretion disks close to the central stars. The exact launching conditions are 
uncertain. We show that a general relation exists between the poloidal and 
toroidal velocity components of a magnetocentrifugal wind at large distances and 
the rotation rate of the launching surface, independent of the uncertain launching 
conditions. We discuss the physical basis of this relation and verify it using a set 
of numerically-determined large-scale wind solutions. Both velocity components 
are in principle measurable from spatially resolved spectra, as has been done for 
the extended low-velocity component (LVC) of the DG Tau wind by Bacciotti et 
al. For this particular source, we infer that the spatially resolved LVC originates 
from a region on the disk extending from ~ 0.3 to ~ 4.0 AU from the star, which 
is consistent with, and a refinement over, the previous rough estimate of Bacciotti 
et al. 

Subject headings: ISM: jets and outflows — magnetohydrodynamics — stars: 
formation — stars: individual: DG Tau — stars: pre-main sequence 
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1. INTRODUCTION 

The magnetocentrifugal mechanism (Blandford & Payne 1982) is widely considered the 
leading candidate for producing jets and winds in star formation (see reviews by Konigl & 
Pudritz 2000; Shu et al. 2000); alternative mechanisms involving thermal and/or radiation 
pressure have been shown to be inadequate. An outflow can in principle be launched by 
torsional Alfven waves, as shown numerically by Shibata & Uchida (1985) and others (see 
Kato, Kudoh, & Shibata 2002 and references therein). These simulations are typically of 
relatively short duration, however, and it is not clear whether a large-scale, sustained outflow 
can be produced by this mechanism alone. 

The magnetocentrifugal model has two popular versions: X-winds (Shu et al. 2000) and 
disk- winds (Konigl & Pudritz 2000). The X-winds are launched from a narrow region close 
to the truncation radius Rx of the Keplerian disk by a stellar magnetosphere. The disk- 
winds are, on the other hand, envisioned to come from a wider range of disk radii outside 
Rx- These two possibilities are not immediately distinguishable observationally, because the 
launching surfaces in both cases are thought to be small, and only the flow regions at much 
larger distances are directly observable at the present. 

In this Letter, we derive a method of inferring the rotation rate (and thus the radius) 
of the wind-launching region on the Keplerian disk from measurable quantities at large 
distances (§2). In §3, we first test the method on a set of numerically-determined large- 
scale magnetocentrifugal wind solutions, and then apply it to the well-studied wind of the T 
Tauri star DG Tau, using data obtained by Bacciotti et al. (2002) with the Space Telescope 
Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST). We discuss our results 
and conclude in § 4. 

2. A METHOD FOR LOCATING THE WIND-LAUNCHING REGION 

2.1. Physical Basis of the Method 

The basic principle of magnetocentrifugal wind launching is well understood. It involves 
open field lines firmly anchored on a rapidly rotating disk, and centrifugal acceleration of 
fluid parcels along the field lines into a high-speed wind. The wind trails behind the disk in 
rotation, generating a toroidal component of magnetic fleld, which exerts a braking torque 
on the disk. It is this magnetic torque that is responsible for extracting both energy and 
angular momentum from the disk and for powering the wind. Since the energy extracted is 
simply the work done by the rotating disk against the magnetic torque (e.g. Spruit 1996), 
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the rate of energy extraction is directly proportional to the rate of angular momentum 
extraction, with the proportionality constant being the angular speed Qq of the disk rotation. 
The extracted energy and angular momentum are initially stored in an electromagnetic 
form. They arc gradually converted into a kinetic form as the flow accelerates. At large 
observable distances, the conversion is nearly complete, and the wind becomes kinetically 
dominated. Along any given field line, the kinetic energy will then be proportional to the 
fluid angular momentum, with the same proportionality constant JIq because of energy and 
angular momentum conservation. 

At any observable location (say of a distance Woo from the rotation axis, where the 
subscript denotes a quantity far from the launching region), the specific kinetic energy and 
angular momentum of the wind are given by the poloidal and toroidal components of the 
fluid velocity, Vp^^o and v^^^o, through ('y|,oo + ^|,oo)/2 and v^^^w^o respectively. Both velocity 
components can in principle be derived empirically from spatially resolved spectra and proper 
motion observations, as has been done for DG Tau (Bacciotti et al. 2002). Once measured, 
they can be used to deduce the rate of disk rotation at the foot point of the field line passing 
through that location, through 

no ^ (1) 

where v^^^ <^ v^^oo is assumed, as is generally true at large distances. The above relation 
provides a simple way to infer the rotation rate of the wind-launching region on the Keplerian 
disk and, if the central stellar mass is known independently, the distance from the star. In the 
next subsection, we will use the steady magnetohydrodynamic (MHD) wind theory to derive 
a refined version of this relation, taking into account of the energy and angular momentum 
associated with fluid rotation at the base of the wind. 



2.2. General Derivation 

It is well known that in a steady, axisymmetric MHD wind, several quantities are con- 
served along a given field hue (Mestel 1968). These include the total specific energy and 
angular momentum 

and the quantity Q = {vtj, — B^Vp/Bp)/zu, which can be interpreted as the angular speed 
at the base of the wind, Qq, where the poloidal flow speed Vp is negligible compared to the 
speed of Keplerian rotation. Here, v denotes velocity, B magnetic fleld, h speciflc enthalpy, p 
mass density, and the gravitational potential. The subscript p refers to a quantity in the 
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poloidal {zu, z) plane of a cylindrical coordinate system {zu, 0, z). We will ignore the enthalpy 
term in the expression for specific energy since magnetocentrifugal winds are dynamically 
cold in general. 

Note from equation (2) that the magnetic contributions to the energy and angular 
momentum are proportional to each other. Since neither of them can be measured directly, 
we get rid of both by constructing a combined quantity 

J = E-VtL = — + ^g- Vtwv^, (3) 

which is also conserved along a field line (see also Lovelace et al. 1986). At the launching 
surface, the wind corotates with the disk at the local Keplerian speed vk,^-, which yields 
J = —TiV^ JI. As the distance increases, the gravitational potential $g decreases quickly, 
and we have approximately 

7;2 _|_ „.2 Q 

v.co ' ^i>,oo r\ '-'2 i w 
2 ^^0^^00-^0,00 ~ -2^^,0- (4) 

Since the Keplerian speed of disk rotation is related to the angular speed through vk.'^ = 
{GM^QqY/^ around a star of mass M*, we finally have 

^00^0,00^^0 - ^(GM.)V3Qf _ <oo + ^loo ^ 

which is the desired equation for determining the disk rotation rate Qq in the wind-launching 
region from measurable quantities at large distances. 

Mathematically, we can define ^ = O^^, and cast equation (5) into a cubic equation 

e - -ao = 0, (6) 

with the coefficients 

02 = , and ao = — . (7) 

'^oo'^</),oo ^'^oo'^4>,oo 

If we let q — —al/Q, r — ao/2 + al/27, and D = + r'^ = a^/A + a2ao/27, then the nature of 
the solution is determined by the value of D. For D > 0, which is always true in our case, 
there are one real and two complex conjugate roots. The real root is given analytically by 

e = (r + ^f'^ + (r - ^^Df'^ + 02/8. (8) 



Once ^ is determined, one can obtain the angular speed through Qq — and infer the 
wind-launching radius through — (GM^/Qq)^''^- 
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Astrophysically interesting magnetocentrifugal winds are probably "fast" in the sense 
that they have enough energy to chmb out the potential well of the central star easily (instead 
of barely). In such a case, we expect the kinetic energy of the wind (the third term of eq. [5]) 
to be substantially greater than the gravitational binding energy at the launching surface 
(which is essentially the second term of eq. [5]). Noting again that typically v^^oo ^ "^p.oo, 
we recover from equation (5) the simple relation (1), which can be cast into a more practical 
form 

for locating the wind-launching region on the disk. 



3. METHOD VERIFICATION AND APPLICATION 

3.1. Method Verification using Numerical Solutions 

We have tested the method of inferring the rotation rate of the wind-launching surface 
on a set of numerically-determined magnetocentrifugal wind solutions (J. Anderson et al. 
2003, in preparation). These solutions are obtained using the Zeus3D MHD code (Clarke, 
Norman, & Fiedler 1994) assuming axisymmetry. The simulation setup and code modifi- 
cations are described in Krasnopolsky, Li, & Blandford (1999). Basically, we prescribe on 
the Keplerian disk around a solar-mass star a distribution of open magnetic field, and load 
onto the field a mass fiux at a low speed (typically 10% of Keplerian). The slowly- moving 
wind material is accelerated centrifugally along the rotating field lines into a high-speed fiow. 
The wind gradually sweeps the ambient medium out of the simulation box, chosen to have 
a 100 X 100 AU size, and settles into a steady state. The properties of the steady wind at 
large distances are determined uniquely by the conditions on the launching surface. They 
are used for our method verification. 

The results of applying our method to two representative simulations are presented 
in Table 1. The simulations differ only in the mass loading rate of the wind; model L 
corresponds to a "light" wind, with My^ = lO~^M0yr~^, typical of T Tauri winds, and 
model H corresponds to a "heavy" wind, with = lO^^Moyr^^, which is more extreme 
(Edwards, Ray, & Mundt 1993). For each model, we first read out from the solution the 
two velocity components Vp^^ and v<^_oo at four selected, observable locations with known 
cylindrical radius zuoo, and then compute the corresponding angular speed at the launching 
surface from both the simple relation (1), ^o,a, and the refined relation (5), Qo,b- Comparing 
the predicted values with the exact Qq, we find that the simple relation (1) works well for 
the light wind, but under predicts the rotation rate by a factor of two for the more extreme. 
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heavy wind. The large discrepancy in the latter case comes from the neglect of the fluid 
energy and angular momentum at the base of the wind in deriving relation (1), which are 
comparable to their magnetic counterparts for the heavy wind. The fluid contributions 
are included in relation (5), which yields a rotation rate accurate to within about 5% in 
both cases. The small errors can be further reduced if the (slow) speed with which the 
wind solution is initiated and the gravitational binding energy at the observing location are 
accounted for properly, which can be done. Such refinements are not warranted, however, 
given the large uncertainties in real observational data, to which we now turn. 



3.2. The Case of DG Tau 

Observational data needed to locate the wind-launching region on the disk arc available 
for DG Tau, thanks to the spatially resolved spectra of several forbidden lines obtained by 
Bacciotti et al. (2000) using HST/STIS. The STIS slit was placed along the axis of the fiow 
and at three pairs of locations at distances 10, 20, and 30 AU from the axis (one on each side). 
Two distinct velocity components are detected, with the high-velocity component (HVC) 
concentrating near the axis and the low-velocity component (LVC) being more laterally 
extended. Bacciotti et al. (2002) analyzed the extended LVC in detail, and derived the 
radial (line-of-sight) velocity ^^,00 at four positions along each slit, labeled 1 through IV in 
their Fig. 7. They also found tentative evidence for rotation from the difference in radial 
velocity between the regions displaced symmetrically with respect to the axis. The axis is 
inclined at ^ ~ 38° to our line of sight (Eisloffel & Mundt 1998), which allows us to estimate 
the poloidal speed Vp^oo — Vr^oo/ cos9, assuming a fiow predominantly parallel to the axis. 
This assumption is expected to break down badly in region I, which is the closest to the disk; 
it should be accurate to ~ 50% or better in the other regions. The estimated values of Vp^^ 
in regions II through IV are listed in Table 2, along with the deprojected toroidal velocities 
computed using the numbers given in Table 1 of Bacciotti ct al. (2002), which are uncertain. 
The angular speed Qq in Table 2 is derived from equation (5), and the wind- launching radius 
zuq is computed for a stellar mass M^, = 0.67 (Hartigan et al. 1995). 

It is clear from Table 2 that the LVC of the DG Tau wind comes from a range in 
disk radius, from ~ 0.3 to ~ 4AU. The spread in launching radius can be understood 
qualitatively as follows: in each of the three regions, the flow located closest to the axis (at 
the 10 AU distance) has the highest poloidal speed (and thus the highest speciflc energy) 
and the lowest, or close to the lowest, toroidal speed (and thus the lowest speciflc angular 
momentum given its small distance). To extract the highest energy against the weakest 
torque on the disk (corresponding to the lowest angular momentum), the field line passing 
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through the innermost location in each region must be driven at the fastest rate at the 
foot point, i.e., must be anchored closest to the central star. Detailed calculations confirm 
this expectation, as shown in Fig. fig:dgtau, where straight lines connecting the observing 
locations and their originating points on the disk are plotted. These lines give a crude 
indication of the actual streamlines. They are well behaved, with two exceptions: the lines 
associated with the two outer locations in region IV cross other lines, which is problematic. 
However, region IV appears to be perturbed by a limb-brightened bubble-like structure 
(Bacciotti et al. 2000), which may have introduced additional asymmetry to the flow with 
respect to the axis, and compromised the estimates of toroidal speed and thus the foot point 
locations. 

We reiterate that our values of Vp^^o are estimated from the radial velocities at the peak 
of forbidden lines, which sample only the component of poloidal velocity parallel to the axis, 
Vz,oo, in the likely case of peak emission originating from the tangent point of our line of 
sight with the wind. Since the component is unmeasured, our values of Vp^^ are actually 
underestimates. This effect will be most pronounced for the least collimated part of the flow, 
where Vp^oo can be underestimated by as much as ~ 50%. Taking this effect into account 
could move the outermost launching point inward by a factor up to two. 

With the rotation rate Qo determined, we can estimate the Alfven radius vja by equating 
the fluid angular momentum at large distances, 'cUooV(p,oo, with the total angular momentum, 
L = Qow\ (e.g. Spruit 1996). The estimates for wa are listed in Table 2 and the correspond- 
ing "Alfven points" are plotted in the right panel of Fig. flg:dgtau. The Alfven radius turns 
out to be ~ 1.8 to ~ 2.6 times the foot point radius (except for the "streamline" passing 
through the bubble-affected outermost location in region IV; see Table 2). The typical ratio 
of wa/'^o ~ 2 implies that roughly 1/4 of the material accreted through the disk is ejected 
in the extended LVC of the DG Tau wind as per the relationship M^j/M^ ~ (^o/'!^a)^ 
(Pelletier & Pudritz 1992). 

4. CONCLUSIONS AND DISCUSSION 

We have developed a method of locating the launching region of T Tauri winds, if such 
winds are driven magnetocentrifugally, as is widely believed. Our method relies on the fact 
that the energy and angular momentum in the wind are extracted mostly by magnetic fields 
from the rotating disk, and they are related by the rate of disk rotation because the energy 
extracted is the work done by the rotating disk against the magnetic torque responsible for 
the angular momentum extraction. Since most of the wind energy and angular momentum 
at large, observable distances are in the measurable kinetic form, they can be used to infer 
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the disk rotation rate in the wind launching region. Applying this method to the wind of 
DG Tau, wc find that its spatially resolved LVC comes from a region on the disk extending 
from ~ 0.3 and ~ 4AU from the central star. This range brackets the rough estimate of 
~ 1.8 AU by Bacciotti et al. (2002). It strengthens the notion that the LVCs of T Tauri 
winds are driven from relatively large disk radii (Kwan & Tademaru 1988). 

The DG Tau wind has a HVC with a radial velocity of ~ 220 km s"^ (Pyo et al. 2003). 
Where this component is launched is less certain. It is unresolved in the transverse direction 
in the HST observations (Bacciotti et al. 2000), implying a half-width < 5AU. Unless 
rotating at a speed much greater than that inferred for the LVC (~ lOkms"-"^), the HVC 
must be driven magnctoccntrifugally from a disk region with radius on the order of 0.1 AU 
or smaller, according to equation (9). The small launching radius is indicative of an X-wind 
origin (Shu et al. 2000). To draw a more quantitative conclusion on the launching radius 
of the HVC, its rotational speed must be measured, which may become possible with the 
advent of optical interferometry using large, ground-based telescopes. 

Support for this work was provided in part by NASA grants NAG 5-7007, 5-9180, 5- 
12102 and NSF grant AST 00-93091. 
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Table 1. Predicted f^o for Model Calculation 
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Cylindrical radius (AU) 



Fig. 1. — Calculated "streamlines" for DG Tau. The left panel shows the observation points 
from Region II (sohd hnes), III (dashed), and IV (dash-dotted) of Bacciotti et al. (2002) 
connected to the calculated foot points of the flow. The right panel is a blow-up of the inner 
region of the flow, showing where the flow originates from the disk. Also shown on the right 
panel is the location of the Alfven surface along each "streamhne" as filled triangles. 
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Table 2. Predicted Qq, vjq and wa for LVC of DG Tau Wind 
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